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ABSTRACT 

This is a review of our experimental and theoretical studies relating to equilibrium and stability of 
liquid drops, typically of low viscosity, levitated in air by a sound fields, 3. The major emphasis here is 
on the physical principles and understanding behind the stability of levitated drops. A comparison 
with experimental data is also given, along with some fascinating pictures from high-speed 
photography. One of the aspects we shall deal with is how a drop can suddenly burst in an intense 
sound field; a phenomenon which can find applications in atomization technology. Also, we are 
currently investigating the phenomenon of suppression of coalescence between drops levitated in 
intense acoustic fields. 


INTRODUCTION 


In the past two decades, or so, acoustic levitation has gained prominence as a way to 
manipulate liquid samples in the microgravity environment of Space [4,5, 6, 7], This is where the 
interest in the study ot interactions between drops and sound arises. The central question is what 
happens when a liquid drop is introduced into a standing sound field. 


The simplest interaction between the drop and the sound field is the force on the 
drop. If the drop surface is assumed to be rigid, then the force on the drop [8] is the summation of the 
normal force acting on the surface due to acoustic radiation pressure <8p>, where [9,10] 
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p ' is the acoustic pressure, u is the particle velocity, p 0 and co are the density and speed of sound of 
the ambient air, and <> denotes average over an acoustic cycle. 


A more complex kind of drop-sound interaction is the deformation of the liquid drop due to the 
acoustic radiation stress. This is resisted by the drop through its surface tension. The static 
deformation is described by the Young-Laplace equation 


P 1 -(8p) = oV.n , (2) 


where o is surface tension, Pi is the internal liquid pressure, and V n represents the total curvature of 
the drop surface [11]. Let us consider a drop, of spherical radius R s> to be at the pressure node of a 
plane wave p '=A sinkz coscot, where A is the pressure amplitude, k is the wavenumber, co is the 
angular frequency, z is the axial position and t is the time. For small amplitude A, the shape of the drop 
is given by [12,13] 
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where B a =R s A 2 /ap 0 c 0 2 is the acoustic Bond number, a=kR s <l represents the spherical radius of the 
drop scaled by the acoustic wavelength A^2rc/k, and P 2 is a Legendre polynomial. Moderately larger 
static deformations have been studied numerically by expansion of the drop shape in Legendre 
polynomials [14]. Also, acoustic deformation can be an excitation source for drop oscillations, if the 
deforming sound field is modulated in amplitude, near a resonant frequency of the capillary oscillations 
of the drop (12,15). 


559 



EXPERIMENTAL OBSERVATIONS 


A common approach to ground-based study of liquid drops is to levitate a millimeter-sized drop 
in a single-axis acoustic levitator, which consists of a transducer at the bottom and a reflector at the 
top, providing a standing sound field that vibrates vertically and falls off laterally. It is known that the 
levitated drop ends up disintegrating in a burst if the sound pressure level (SPL) is raised too high 
[16], Our goal is to understand why, when, and how the burst can occur. The first step is to study 
the phenomenon in detail, using high-speed photography. The observations are described as follows 
[ 1 , 2 ]. 


All drops levitated in a sound field become more flattened with increasing SPL, progressively 
attaining oblate spheroidal, and then flat disk shapes. Their subsequent evolution and stability can be 
differentiated based on their initial size. Two limiting scenarios can be identified. Whenever the drop is 

acoustically small (a<0.3), flattening cannot exceed a critical point, beyond which the drop 
disintegrates. The manner of disintegration is a sudden horizontal expansion, in which the drop is 

pulled apart at the periphery (Fig. 1). Acoustically large drops(a>0.5), on the other hand, can sustain 
equilibrium beyond the critical point by actively lowering the sound pressure level in the levitator. A 
large flattened drop is an efficient scatterer of the acoustic wave and can induce a sufficiently large 
resonance frequency shift in the system, causing the levitator to lose power [17]. The SPL in the 
levitator first increases with flattening, as expected, but then decreases with flattening beyond the 
critical point. Here, the drop becomes concave and eventually drastically flattened to the point where 
the drop has a meridional cross-section resembling a dogbone, with a membrane in the middle 

surrounded by a donut-shaped periphery (fig. 2). Drops of intermediate size (0.3< a<0.5) lose 
equilibrium beyond the critical point and show sudden horizontal expansion, but they return to 
equilibrium following sufficient lowering of the sound pressure level through resonance frequency shift. 

On continued flattening, the central membrane of a large drop may eventually buckle upward, 
with the ring-shaped periphery contracting and closing, leading to violent shattering (Fig.2 f,g,h). The 
buckling is always upward, such that it has to be related to the presence of gravity. For all low- 
viscosity drops, like water drops, as the central membrane vibrates in the sound field capillary ripples 
are excited parametrically on it. Here, unless buckling instability occurs first, the ripples become so 
intense that satellite droplets are ejected vertically in both directions. This phenomenon is 
conventionally referred to as atomization. For small drops, on the other hand, atomization takes place 
near the edge, during the sudden horizontal expansion, and is inherently nonsteady. A highly viscous 
drop does not atomize, by this route, because any ripple would be heavily suppressed by viscosity. 


STATIC EQUILIBRIUM 

To understand these complex drop-sound interactions, the natural thing to do is to understand 
the equilibrium of the system, then find out when the system becomes unstable [1 ,2,3]. The poles of 
the drop are stagnation points of the oscillating flow of the sound field. In these regions, according to 
eq.(1), the radiation stress comes mostly from the acoustic pressure and is positive. At the equator, on 
the other hand, the particle velocity is most intense. The radiation stress here comes mostly from the 
particle velocity and is negative. The drop shape is determined by the compressive stress at the 
poles, the suction stress at the equator, together with surface tension, and the constant internal 
pressure inside the drop. When the central membrane is thin enough, it vibrates with sound a little. In 
our cases of interest, the vibration is not strong enouah to significantly violate the rigid wall boundary 
condition. The membrane vibration leads to a Bernoulli correction in the internal pressure. The 
equilibrium shape, as a consequence of all these forces, changes from an oblate spheroid to one with 
flat poles, and then to one with concave poles. 

Mathematically, the difficulty with the problem of a drastically deformed drop is that the drop 
shape depends on the wave, but the wave also depends on the shape. To handle the problem in a 
self-consistent manner, we have used a boundary integral method for solving the Helmoltz equation for 
the wave [3,18], given the shape. Then we use the resulting radiation pressure to solve for the shape 
using eq.(2). The two procedures are repeated alternatively until both the wave and the shape 
converge to their respectively forms. 
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For a given drop size a, if we plot the acoustic Bond number B a versus the non-dimensional 
equatorial radius R* (R*=R/R S , R: equatorial radius), we find that B a rises with R* to a maximum, and 
then decreases with R*, schematically as shown by the solid curve in Fig.3. Correspondingly, the 
SPL increases to a maximum and then decreases with flattening, in agreement with the observation for 
large drops. A large drop can bring about a downturn in SPL because it can easily induce the 
necessary resonant frequency shift. Therefore it stays on the solid equilibrium curve until buckling 
instability, or atomization, occurs. An intermediate sized drop has to expand sufficiently, before it can 
bring about the necessary lowering in sound intensity, to be on the equilibrium curve. Therefore the 
drop follows curve 1 (fig. 3) which jumps from the equilibrium curve, but returns to it later. A small drop 
on the other hand, cannot bring about a significant power loss by resonant frequency shift during 
expansion. Hence it follows a curve like 2 (fig. 3) to disintegration. 

A typical B a -R* curve for a large drop is shown in Fig.4, showing good agreement with 
experiment. When the central membrane of a large drop vibrates, the internal pressure in the 
membrane cannot be lowered much because the surface curvature of a thin membrane is close to zero. 
Thereby, the internal pressure has to be close to the external radiation pressure. The Bernoulli effect, 
due to membrane vibration, is manifested in the relative raising of the liquid pressure in the ring-shaped 
periphery. This tips the force balance at the edge, such that the drop loses equilibrium at a lower 
flattening than that for which central membrane does not vibrate. If the drop is viscous, vibration of the 
membrane is spread into the periphery too, such that Bernoulli effect is effectively nullified and can be 
ignored. The point here is that a viscous drop can be flattened more than a inviscid drop because of 
the absence of Bernoulli effect due to membrane vibration. This is consistent with experimental 
observations. 

The actual occurrence of the downturn of the B a -R* cun/e depends on the active involvement of 
the levitator. In practice, we may be more concerned with drops that are too small to cause any power 
loss, or a sound system that is too strong to allow for a significant power loss. Hence the peak of a B a - 
R* curve is important because it represents the maximum sound intensity and maximum deformation a 
drop of size a can sustain. In Fig.5, we plot the maximum B a , denoted by B acr , and the corresponding 
R*, denoted by R* cr , versus a. Experimental data are also displayed for comparison, showing good 
agreement between the two. 


INSTABILITIES 

On the downturn side of the equilibrium curve (fig. 3), a large drop will lose equilibrium at the 
edge, beyond a certain degree of flattening. But before the drop reaches this point, instabilities will 
occur. These can be either: (i) the buckling instability [1] or, (ii) the parametric instability that gives rise 
to the ripples and the subsequent atomization [2, 19]. 

Buckling instability is explained as follows. When a drop is drastically flattened such that the 
membrane is thin enough to be forced into vibration, like a drumhead by the sound field, a dynamic 
pressure arises. The latter is given by p<v2>/2, where p is the density of the liquid, and v is the 
vibration velocity of the membrane. This pressure tends to push the liquid toward the periphery. If 
equilibrium is to be maintained, the hydrostatic pressure in the membrane has to become lower, as 
dictated by Bernoulli law. But the membrane surfaces are almost flat, such that there is limited room for 
readjustment. The membrane becomes more susceptible to acoustic vibration as it gets thinner with 
flattening. Ultimately, there is a point beyond which the hydrostatic pressure in the membrane can no 
longer withhold the dynamic pressure. If the membrane is perfectly symmetric about its mid-surface, 
there is still the relatively higher hydrostatic pressure in the thick periphery to contain the membrane 
liquid. But in the presence of gravity, the membrane is slightly buckled upward, to start with. In this 
case, the dynamic pressure can cause the membrane to expand and buckle upward, as shown in Fig.2 
(f,g,h). 


For a low viscosity drop, capillary ripples can be excited on the membrane through parametric 
instability, when the membrane vibrates in the sound field (frequency to). In general, there can be two 
wave modes: a symmetric wave (frequency oo s ) in which the two surfaces vibrate in opposite 
directions, and an antisymmetric wave (frequency co a ) in which they vibrate in the same direction [20]. 
For a vibrating membrane, of a given thickness, the wavelength of the disturbance is determined as a 
consequence of parametric instability [2,19]. Parametric excitation of a wavelength can occur only when 
the wavelength satisfies the synchronism condition: co=(o s +co a . When the membrane is thick 
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co a =co s =o a/2, and the waves are known as Faraday ripples. When the ripples grow violent, they spit 
out satellite droplets as in conventional atomization. When the membrane is thin, co a is much larger than 
o) s . A thin film will shatter as ripples grow because it will snap wherever its two surfaces come into 
contact. If a drop is very viscous, then ripples are suppressed. The easier route for a viscous drop to 
fragment is through a sudden horizontal expansion, in which the membrane is stretched so thin that it 
has to rupture and shatter due to van der Waals forces. 

CLOSING COMMENTS 

In this article, we have addressed the issue of equilibrium and stability of a liquid drop, when 
introduced into an intense sound field. This complex problem has been primarily presented in a 
physical way, substantiating it with clear photographic evidence and a simple theoretical framework 
employing the acoustic Bond number and the drop equatorial radius (B a -R* curve). Presently, we are 
trying to understand the mechanism of suppression of coalescence of drops in intense acoustic fields. 
Normally, when two drops are introduced into a levitating field , they coalesce to form a single drop. 
However, whenever the intensity of the acoustic field is high, the two drops do not coalesce, but 
instead form a stable doublet (fig. 6 a,b). They seem to contact along a line, and their contacting 
surfaces display a knife-edge feature. The reason for this stability has probably to do with the nature 
of the vibrating air flow in the narrow strait between the drops, and needs to be well understoodsi . The 
stable doublet can be rotated as a rigid body, without the two drops separating. The drops coalesce 
only when the intensity of the levitating field is reduced, thereby making the drops less flattened and 
more spherical. 
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Fig. 1 The side view of a small water drop (a=0.2): (a) 
being flattened, (b) its rim getting sharper, (c) rim 
developing into a knife-edge, (d,e) expanding 
horizontally and instability occurring along the 
outermost edge, and (f) final breakup. 


Fig. 2. The evolution of a large water drop (a=0.8) 
levitated in an acoustic field with increasing intensity: 
(a) being flattened and slightly dimpled, (b,c) 
becoming concave, (d,e) ripples formed on the 
central membrane, spreading throughout the latter, 

(f) buckling upward, (g) ballooning with two jets 
shooting in opposite directions at the bottom part that 
has just been closed, and (h) final breakup. / 



Fig. 3. Responses of a levitated drop to a continuous increase in input sound intensity. 
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Fig. 4. Plot of Ba vs R* for a water drop with a=0.74. 
The crosses are experimental data and the solid 
squares are theoretical data. The dotted curves are 
asymptotic forms given by Ba=Bm*/R*, where 
Bm*=3.52 for the theory and 3.95 for the experiment. 


Fig. 5. The maximum Ba of the Ba-R* curve, denoted 
by Ba.cr , and the corresponding R*, denoted by R*cr 
, for a range of a. Theoretical data are given by 
squares and experimental data are given by crosses. 



Fig. 6 . Side (a) and top (b) views of a stable non-coalescing doublet in intense acoustic field. 
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